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RAYLEIGH-SOMMERFELD  DIFFRACTION  YS  FRESNEL-KIRCHHOFF, 
FOURIER  PROPAGATION,  AND  POISSON’S  SPOT 


1.  INTRODUCTION 

The  boundary  conditions  imposed  on  the  diffraction  problem  in  order  to  obtain  the  Fresnel-Kirchhoff 
(FK)  solution  are  well-known  to  be  mathematically  inconsistent  and  to  be  violated  by  the  solution  when 
the  observation  point  is  close  to  the  diffracting  screen  [1-3].  These  problems  are  absent  in  the  Rayleigh- 
Sommerfeld  (RS)  solution.  The  difference  between  RS  and  FK  is  in  the  inclination  factor  and  is  usually 
immaterial  because  the  inclination  factor  is  approximated  by  unity.  But  when  this  approximation  is  not  valid, 
FK  can  lead  to  unacceptable  answers.  Calculating  the  on-axis  intensity  of  Poisson’s  spot  provides  a  critical 
test,  a  test  passed  by  RS  and  failed  by  FK.  FK  fails  because  (a)  convergence  of  the  integral  depends  on  how 
it  is  evaluated  and  (b)  when  the  convergence  problem  is  fixed,  the  predicted  amplitude  at  points  near  the 
obscuring  disk  is  not  consistent  with  the  assumed  boundary  conditions. 

Poisson’s  spot,  also  known  as  the  spot  of  Arago,  is  the  name  given  to  the  bright  on-axis  spot  behind  a 
circular  obscuration  illuminated  by  a  plane  wave:  on  the  axis  of  the  disk,  all  light  diffracted  at  the  rim  of  the 
disk  arrives  in  phase  and  interferes  constructively.  Consequently,  even  for  angles  approaching  90°,  i.e.,  for 
observation  points  close  to  the  disk,  diffraction  can  result  in  a  significant  intensity.  (In  the  overwhelming 
majority  of  optics  problems,  diffraction  at  angles  of  more  than  a  few  degrees  leads  to  vanishingly  small 
intensities  and  can  be  safely  ignored.)  Treatments  of  Poisson’s  spot  that  make  the  Fresnel  approximation  [4,5] 
do  not  apply  to  the  region  close  to  the  disk.  Treatments  that  do  apply  close  to  the  disk  use  the  RS  integral 
[6]  or  its  predecessor,  the  Rayleigh  integral  [7].  The  RS  integral  is  derived  from  the  Rayleigh  integral  with 
the  additional  assumption  that  the  wavelength  of  the  light  is  small  compared  to  the  geometric  dimensions  of 
the  problem,  an  assumption  that  is  made  throughout  this  report.  I  refer  to  a  calculation  as  exact  to  indicate 
that,  while  the  short-wavelength  approximation  has  been  made,  the  Fresnel  approximation  has  not.  The  other 
fundamental  approximation  made  here  is  the  scalar  wave  approximation,  which  means  that  the  results  apply 
better  to  acoustics  than  to  optics,  a  point  that  will  be  reconsidered  in  the  concluding  remarks  of  Section  5. 

Fourier  propagation  provides  an  alternate  means  of  handling  diffraction  problems.  It  can  be  used  in 
computer  software  to  calculate  beam  propagation  problems,  but  in  this  report,  Fourier  propagation  theory 
is  used  to  derive  solutions  in  integral  form  for  the  2-D  (long  slits  and  strips)  and  3-D  (arbitrarily  shaped 
apertures)  diffraction  problems.  Standard  Fourier  propagation  reproduces  the  RS  diffraction  integral. 

2.  THE  BASIC  DIFFRACTION  PROBFEM 

Solving  the  basic  diffraction  problem  requires  finding  a  solution  to  the  Helmholtz  equation  for  a 
propagating  wave  encountering  a  partially  obscuring  planar  screen.  The  Helmholtz  equation  is 

[v2  +k2^U(x,y,z)  =  0,  (1) 
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where  k  =  2k  fk  and  U  describes  the  amplitude  and  phase  of  the  wave.  U  is  a  scalar,  so  only  scalar  diffraction 
theory  is  addressed.  The  boundary  condition  imposed  on  the  solution  to  this  differential  equation  is  the  effect 
of  a  diffracting  screen  in  the  z  =  0  plane.  Denoting  by  T  the  parts  of  the  screen  that  are  transmissive  and  by 
B  the  parts  that  block  the  beam,  the  boundary  conditions  used  for  the  RS  and  FK  solutions  are 


RS  and  FK: 


FK  only: 


U(x,y,0)  =  U0(x,y,0) 
U(x,y,  0)  =  0 


dU(x,y,z) 


dz 

dU(x,y,z ) 


dU0(x,y,z) 


z= 0 


dz 


z= 0 


dz 


=  0 


z=0 


for  (x,y)  e  T , 
for  (x,  y)e  B, 

for  (x,y)eT, 
for  (x, y)&B, 


(2) 


where  Uq  describes  the  incident  wave  and  dU/dz  is  the  derivative  of  U  normal  to  the  diffracting  plane. 

With  reference  to  Fig.  1,  the  RS  diffraction  integral  [1,2]  for  U  at  distance  z  behind  an  aperture  in  a 
planar  mask  is 


URs(,x',y',z)  =  “  J  U0(x,y,0)CXp  ^  cos  yylxdy ,  (3) 

^  Area  r 

ry  9  9V. 

where  r  =  [(jc'  -  x)  +  (y'  -y)  +  z  ]  is  the  distance  between  (x,y,  0)  and  (x',y\ z),  y  is  the  diffraction  angle 
at  point  (x,  y,  0),  i.e.,  the  angle  the  diffracted  ray  makes  with  the  normal  to  the  plane  (not  with  the  direction 
of  the  incoming  wave)  and  the  integral  is  over  the  area  of  the  aperture.  The  FK  integral  [1,3]  is 

UFK(x',y',z)  =  ~Y  J  £/0(x,y,0)eXP  ^  (cos  ^  +  cos  y)dxdy,  (4) 

^  Area  r  ^ 

where  'C  is  the  incidence  angle  at  point  (x,  y,  0).  The  exp (ikr)/r  factors  in  Eqs.  (3)  and  (4)  express  Fluygens’ 
principle:  each  point  in  the  aperture  acts  as  a  source  of  spherical  waves  that  combine  to  give  the  diffraction 
pattern.  The  cosine  factors  are  called  the  inclination  factors  and  constitute  the  only  difference  between  RS 
and  FK.  For  most  diffraction  problems,  the  inclination  factor  is  approximated  by  unity,  causing  the  difference 
between  RS  and  FK  to  disappear. 

Goodman  [1]  gives  a  succinct  derivation  of  both  RS  and  FK  and  discusses  the  difference  between  them. 
The  derivations  use  different  Green  functions,  which  require  different  boundary  conditions  to  reduce  the 
Green’s  theorem  integral  to  the  familiar  diffraction  integrals  given  in  Eqs.  (3)  and  (4).  The  basic  problem  with 
FK  is  that  the  Green  theorem  integral  can’t  be  evaluated  unless  the  values  of  both  U  and  dU/dz  are  assumed 


Fig.  1  —  The  basic  diffraction  problem:  a  point  source  p 
illuminates  an  aperture  and  produces  a  diffraction  pattern 
on  an  observation  screen.  The  distance  r  is  the  distance  from 
a  point  in  the  aperture  to  a  point  on  the  screen.  The  angles  £ 
and  x  (not  labeled  in  the  figure)  are  the  angles  the  incident 
and  the  diffracting  rays  make  with  the  z  axis. 
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to  be  zero  on  the  obscuring  part  of  the  diffracting  screen,  but  we  know  from  analytic  function  theory  that  if 
both  U  and  dU/dz  are  zero  over  any  region,  then  (7=0  everywhere.  Thus,  the  FK  solution  cannot  be  fully 
mathematically  consistent  and  must  therefore  be  suspected  of  not  always  (at  least)  giving  the  right  answer. 
In  Section  3, 1  use  Poisson’s  spot  to  show  explicitly  that  the  FK  solution  can  violate  the  boundary  conditions 
that  were  assumed  for  its  derivation. 


3.  CALCULATING  THE  INTENSITY  OF  POISSON’S  SPOT 


For  the  simple  case  of  a  plane  wave  impinging  normally  on  a  circular  disk  of  radius  a  (Fig.  2), 
cost  =  1  and  cosy  =  z/r.  I  calculate  the  on-axis  amplitude  and  intensity  behind  the  disk  by  doing  so  for  an 
annular  aperture  and  letting  the  outer  radius  of  the  annulus  approach  infinity: 


U(0,0,z)  =  - 


p—>°° 

U0  J 

p =a 


exp  (ikr)  f 


q+c2- 
V  r 


27Tpr/p, 


(5) 


where  q  =  0,  c2  =  2  for  RS;  q  =  c2  =  1  for  FK;  and  p  is  the  radial  coordinate  in  the  (x,  y)  plane.  In  order  to 
evaluate  Eq.  (5)  without  using  the  Fresnel  approximation,  I  follow  Sommerfeld  [2]  and  Harvey  et  al.  [6]  in 
changing  the  variable  of  integration  from  p  to  r.  On  the  z  axis,  r=(z+  p  )  ,  so  a  <  p  <  oo  =>  r0  <  r  <  oo, 
where  r0  =  (z2  +  a2)'/2.  Now  r2  =  z2+  p2,  so  rdr  =  pc/p,  and  the  integral  in  Eq.  (5)  can  be  put  in  the  right  form 
to  be  evaluated  via  the  Sommerfeld  lemma  given  in  Appendix  A: 


r—>°°  j  (  7 

U(0,0,z)  =  -HcUq  J  exp  (ikr)—  C|  +  c2  -  I  dr 

2  V  r 

'o  v 


=  -ikU, 


—  exp  (ikr)  — 
ik  2 


r— 

f  1  1 

Ct  +  c2  - 

+o 

72 

l  rj 

r=ro 

\k~  ) 

-Un 


q  +  c2- 


r0 


exp  (ikr0  )-U0^-  exp  [ik(r  — >  °°)] . 


(6) 


The  RS  version  of  Eq.  (6)  is 


URS(0, 0,  z)  =  U0  —  exp  (; ikr0 ) , 
r() 


(7) 


while  the  FK  version  is 


Ufk(0,0,z)  =  Uq 


1 

2 


V  roJ 


exp  ( ikr0)-U0 


1 

2 


exp  [  ik(r  — >  oo)] , 


(8) 


which  shows  that  the  FK  integral  fails  to  converge  in  this  case.  The  reader’s  attention  is  called  to  the  fact  that 
if  a  is  set  to  zero  (i.e.,  r0  =  z),  there  is  no  obscuration  and  the  right  side  of  Eq.  (8)  should  be  just  U(fxp(ikz) 
-  which,  because  of  the  second  term,  it  isn’t!  The  FK  integral,  when  evaluated  in  this  straightforward  way, 
doesn’t  give  an  acceptable  answer. 
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Fig.  2  —  The  diffraction  problem  for  Poisson’s 
spot;  a  plane  wave  falls  on  a  circular  disk.  The 
distance  r  =  (x2  +  y2  +  z2)l/2  is  the  distance 
from  a  point  in  the  plane  of  diffraction  to  a 
point  on  the  axis  of  the  disk. 


Eliminating  the  second  term  in  Eq.  (8)  can  be  done  in  various  ways.  An  artificial  way  would  be  to  make 
the  outer  edge  of  the  annular  aperture  an  ellipse  instead  of  a  circle.  Then  the  rays  diffracted  from  this  edge 
do  not  arrive  on  the  axis  in  phase  and  do  not  interfere  constructively.  But  the  most  sensible  way  is  to  impose 
Babinet’s  principle  as  a  separate  requirement  (separate,  because,  as  we  have  just  seen,  FK  doesn’t  satisfy 
it  unless  the  integral  is  done  the  right  way).  Babinet’s  principle  requires  that  the  sum  of  the  obscuration 
diffraction  pattern  (Fig.  2)  and  the  aperture  diffraction  pattern  (Fig.  2  with  the  transmitting  and  blocking 
regions  reversed)  be  the  uninterrupted  plane  wave:  Uoh  +  Uap  =  U^cxpU/tz).  Thus  Uob  =  (70exp(/kz)  -  Uap , 
which  can  be  seen  from  Eqs.  (5)  and  (6)  to  be  Eq.  (8)  without  the  second  term: 


UFK(0,0,z)  =  U o  exp  (ikz)+~U0j 

A  o 


exp  (ikr  )  1 


1  +  -  2npdp 
V  r , 


1 

=  Un  — 


1+*- 

>0. 


exp  (ikt\) ) . 


(9) 


The  remaining,  and  more  serious,  problem  with  Eq.  (9)  is  that  it  doesn’t  satisfy  the  boundary  condition 
under  which  it  was  derived.  As  z  ->  0,  we  should  find  (7(0,  0,  z)  ->  0  in  Eqs.  (7)  and  (9),  since,  from  Eq.  (2), 
that  is  the  boundary  condition  originally  assumed.  Eq.  (7)  satisfies  this  condition,  while  Eq.  (9)  doesn’t.  Also, 
omitting  the  exp  (ika)  phase  factor,  dU/dz=  Uq/u  at  z  =  0  for  RS  and  half  this  for  FK.  The  non-zero  value  of 
dU/dz  is  consistent  with  RS,  for  which  only  U  need  be  zero  at  the  disk,  but  not  for  FK,  which  began  with 
the  additional  boundary  condition  dU/dz  =  0. 


Eqs.  (7)  and  (9)  are  squared  to  obtain  intensity  and  are  plotted  in  Fig.  3.  The  predicted  intensities  begin  to 
differ  appreciably  at  z/a  «  4,  where  the  diffraction  angle  is  %  ~  15°.  Note  that  the  exactly  on-axis  intensity  of 
Poisson’s  spot  doesn’t  depend  on  wavelength  because  there  is  no  wavelength  dependence  in  the  inclination 
factor.  Wavelength  dependence  enters  in  the  radial  intensity  distribution,  which,  for  the  central  region  defined 
by  rja  «  1 ,  where  r±  is  the  2-D  radius  measured  from  the  z  axis,  can  be  shown  to  be  [4] 


U(rF,z)  =  U(0,0,z)  exp 


A.  2  A 

/7trj_ 

fa 


r  2narF  ) 


(10) 


The  first  zero  of  70  is  at  2nar1lfa  =  2.4,  which  shows  that  the  diameter  of  Poisson’s  spot  is  proportional  to 
A,  so  the  spot’s  area  and  the  power  contained  in  it  are  proportional  to  A2. 

4.  DOING  DIFFRACTION  WITH  FOURIER  PROPAGATION 

In  Section  4.1,  I  invoke  the  basic  principles  of  Fourier  propagation,  and  then  use  these  principles  in 
Section  4.2  to  derive  the  RS  diffraction  integral.  For  ease  of  presentation,  I  address  the  2-D  diffraction 
problem:  long  slits  or  strips  illuminated  by  plane  or  cylindrical  waves  that  can  be  described  by  U(x,  z).  In 
Section  4.1  the  generalization  to  3-D  is  obvious;  for  Section  4.2,  it  is  done  in  Appendix  B. 
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Poisson's  Spot 


Fig.  3  —  Intensity  of  Poisson’s  spot  at  on-axis  distance  z  behind 
a  disk  of  radius  a ,  according  to  the  Rayleigh-Sommerfeld  and 
Fresnel-Kirchhoff  diffraction  theories. 


4.1  Fourier  Propagation 


The  Fourier  transform  over  x  of  U(x,  z )  is  called  its  angular  spectrum,  defined  by 


A(a,z)=  J  U(x,z)  exp  (-ikajc)dx. 

— oo 

Fourier  propagation  rests  on  the  premise  that,  as  shown  by  Goodman  [1]  for  example, 

A(a,z)  =  A(a,0)  exp  lilts]  1-a2. z 


(11) 


(12) 


where  a  is  not  restricted  to  -1  <  a  <  1.  The  Fourier  transform  variable  in  Eq.  (11)  is  afk,  so  the  inverse 
transform  of  A(a,  z)  is 


U(x',z)=  J  A(a,  z)  exp  (ikax')d 


=  J  A(a,0)  exp  liksll-a2z  +  ikax'\  d 


=  J  J  U(x,0)  exp  (-ikax)dx  exp  ( iksl  1  -  a2  z  +  ikax'\  d 


— oo  — oo 


oo  oo 


"a 

J 


=  J  J  t/(x,0)exp 


— oo  — oo 


iksll-a2 . z  +  ika.(x'  —  x) 


dxd 


(13) 


The  second  equality  shows  that  U{x' ,  z)  is  the  sum  of  plane  waves  of  amplitude  A( a,  0),  propagating  at 
angle  0  =  cos* 1  a  with  respect  to  the  x  axis.  When  |a|  >  1 ,  the  plane  waves  are  evanescent,  with  exponentially 
decaying  z-dependence  given  by  exp  [-(or  -  1)  z].  These  waves  do  not  propagate  a  significant  distance  from 
the  aperture  but  are  needed  to  give  a  complete  Fourier  decomposition  of  U(x',  z)  at,  or  near,  z  =  0. 
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Since  the  dependence  on  x'  and  z  is  all  in  the  exponent  on  the  right  side  of  Eq.  (13),  it  is  easy  to  see 

that 


r 


v 


a2 


dx,z 


+ 


U(x',z ) 


U(x',z ) 


-k’U(x',z), 


(14) 


which  shows  that  U(x',  z)  is  a  solution  to  the  2-D  Helmholtz  equation,  and  Eq.  (13)  shows  that  the  solution 
depends  only  on  the  value  of  U  in  the  z  =  0  plane.  For  the  diffraction  problem,  the  standard  procedure  is 
to  use  U()(x,  0)  for  U(x,  0)  in  the  transmissive  parts  of  the  screen  and  zero  in  the  blocking  parts.  (Goodman 
loosely  states  that  “Kirchhoff  boundary  conditions”  are  applied  [1],  but  actually  only  the  RS  conditions  are 
needed.)  Thus,  Fourier  propagation  and  the  RS  integral  are  solutions  to  the  same  differential  equation  with 
the  same  boundary  condition,  hence  must  be  the  same  function,  a  fact  that  is  shown  explicitly  in  the  next 
section. 

4.2  The  RS  Integral  Derived  via  Fourier  Propagation 


Returning  to  Eq.  (13),  interchanging  the  order  of  integration  and  using  the  RS  boundary  condi¬ 
tion, 


U(x,z)  =  | t/0(x,0)j  exp 

x  a 


=  J U0(x, 0)  J  exp  [/foF(a)]  d 
x  a 


dx, 


dx 


(15) 


where  the  last  line  defines  the  function  F(a).  The  a  integral  is  done  by  the  stationary  phase  method:  the 
function  F(a)  is  expanded  in  a  Taylor  series  about  the  point  a0  at  which  its  first  derivative  is  zero.  Using 
the  notation  Fa  =  dF/da,  we  require 


which  can  be  solved  for  a0, 


F<x  (®o ) — 


x  -  x  x  -  x 

«0  =  /  9  9  = 
V(*'-;t)2+r  r2  d 

Therefore 

V1_(X  0  =  — 

r2D 


(16) 


(17) 


(18) 


and 


F(a0)  =  Jl-al  +a0  - — -  =  — 

z  z 


(19) 
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The  second  derivative  of  F(a )  at  a0  is 


^aa(ao)~  I - 

•\/1_ao 


(20) 


Since  the  first-order  term  vanishes  by  construction,  the  Taylor  series  expansion  through  second  order  of 
F(a )  about  a0  is 


1  9 

F(a)  ~  F(a0 )  +  —  Faa  (a0  )(a  -  a0 )  . 


A  standard  Fresnel  integral  is  now  written  in  a  form  that  will  be  useful  below: 


J 


-j-A(u-u0)2 


l-i  1 

V2  VXiA’ 


(21) 


(22) 


where  A  is  positive  definite. 

The  integral  over  a  in  Eq.  (15)  can  now  be  evaluated.  For  ease  of  notation,  drop  the  a0  argument  from 
Faa  and  observe  that  Faa  =  -\Fua\  (Faa  is  always  negative).  Eq.  (21)  is  inserted  into  the  exponent  in  Eq. 
(15)  to  obtain: 


exp  [/foF(a0)]  J  exp 
a 


(a-a0)2 


=  exp  (ikr2D) 
=  exp  (ikr2D) 


1  —  i 


1  —  i 

12 


1 


aa\ 

1  Z 


sjRf20  r2D 


(23) 


Therefore 


U(x',z ) 


J  Uq(x,0)  exp  (ikr2D) 

X 


l-i  1  z 

■n/2  yjh 2D  r2D 


dx 


l-i  r  exp  (ikr7n) 

-j=^  J  U0(x, 0)  )— cos  %2 Ddx- 


(24) 


Eq.  (24)  is  the  2-D  version  of  Eq.  (3)  and  contains  the  2-D  Eluygens’  principle:  each  strip  in  the  aperture  acts 
as  a  source  of  cylindrical  waves,  for  which  amplitude  falls  off  as  The  interested  reader  may  want  to 
make  the  normal  approximations  [cosx2d  “  1  >  r2D  =  z  +  {x'  -  x)2/2z  in  the  exponent,  r2D  »  z  outside  it]  and 
show  that  Eq.  (24)  reduces  to  the  standard  form  that  is  evaluated  with  the  Cornu  spiral. 
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Extending  Eq.  (24)  to  the  3-D  problem  will,  first  of  all,  add  a  dy  to  the  integral.  Keeping  in  mind  that 
all  the  units  of  length  must  cancel  out  on  the  right  side,  inspection  of  Eq.  (24)  suggests  that  the  effect  of  a 
3-D  calculation  is  to  replace  r2Q  and  X21)  by  r  and  x,  to  replace  cylindrical  waves  by  spherical  waves,  and 
to  square  the  factor  outside  the  integral.  This  intuitive  argument  is  confirmed  in  detail  in  Appendix  B,  with 
a  result,  Eq.  (B14),  that  matches  the  RS  form  of  the  diffraction  integral  given  in  Eq.  (3). 

5.  CONCLUSION 

The  fundamental  flaw  in  the  FK  diffraction  integral  and  the  superiority  of  RS  have  been  demonstrated 
with  exact  calculations  of  the  intensity  of  Poisson’s  spot.  Fourier  propagation  has  been  presented  as  an 
alternate  means  of  deriving  the  diffraction  integral.  Compared  to  the  usual  approach  via  Green’s  theorem, 
this  derivation  has  the  advantage  of  rendering  obvious  the  proper  choice  of  boundary  conditions.  It  has  the 
disadvantage  of  requiring  knowledge  of  Fourier  propagation  and,  for  the  3-D  version,  more  complicated 
math,  but  the  2-D  version  is  not  bad  and  the  generalization  to  3-D  by  inspection  is  intuitively  appealing. 

As  noted  in  the  introduction,  the  argument  has  been  confined  to  scalar  wave  theory,  which  will  not  be 
completely  adequate  for  describing  Poisson’s  spot  in  optics  for  points  near  the  disk.  (The  contribution  of 
those  rays  not  polarized  parallel  to  the  diffracting  edge  should  be  multiplied  by  the  sine  of  the  angle  between 
the  polarization  vector  and  the  direction  to  the  observation  point.)  But  it  should  be  entirely  adequate  for  de¬ 
scribing  Poisson’s  spot  in  an  acoustics  experiment  because  acoustic  waves  are  scalar  (pressure)  waves.  Also, 
because  the  wavelength  is  much  longer,  diffraction  phenomena  can  be  more  easily  studied  in  an  acoustics 
than  in  an  optics  laboratory.  This  experiment  was  done  many  years  ago  with  somewhat  equivocal  results 
[8].  With  modem  equipment,  it  should  not  be  particularly  difficult  to  repeat,  and  could  settle  the  conflict 
between  RS  and  FK  diffraction  by  direct  measurement. 
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Appendix  A 


THE  SOMMERFELD  LEMMA 


Following  Sommerfeld  [Al]  and  Harvey  et  al.  [A2],  we  perform  a  series  of  integrations  by  parts  (only 
the  first  two  are  shown)  to  expand  the  integral  of  a  function  multiplied  by  a  complex  exponential  in  a  series 
of  terms: 


b 

\  fix)  exp  (ikx)  dx 
a 


f  O) 

ik 


exp  (ikx) 


f\x) 

- —  exp  (ikx) 

(ik)2 


(ik)1  a 


J  f"(x)  exp  (ikx)dx 


^  ^  exp  (ikx) 
ik 


b 

a 


(Al) 


where  the  approximation  is  justified  if  f(x)  is  a  slowly  varying  function,  or,  equivalently,  in  the  limit  k—>° ° 
(k  — >  0).  When  applied  to  diffraction,  the  approximation  in  Eq.  (Al)  holds  when  the  geometric  dimensions 
of  a  problem  are  large  compared  to  a  wavelength  of  light,  otherwise /(x)  may  not  be  sufficiently  slowly 
varying. 

The  integral  in  Eq.  (Al)  has  exactly  the  form  of  the  ^-frequency  coefficient  of  a  Fourier  series  expan¬ 
sion  of  the  function/(x)  over  the  interval  (a,  b).  This  shows  that  for  large  k,  we  need  know  only  the  values 
of /(a)  and  f(b)  to  find  the  value  of  the  coefficient.  I  have  searched  more  than  a  dozen  Fourier  series  books 
and  have  not  found  this  lemma  in  any  of  them. 
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Appendix  B 


DERIVING  THE  RS  DIFFRACTION  INTEGRAL  VIA  FOURIER  PROPAGATION 

Proceeding  in  analogy  to  Section  4.2,  the  3-D  version  of  Eq.  (15)  is 
U(x',  y',  z ) 


=  J  U0(x,y, 0)  J  exp 

x,y  a,P 


f  r 

ikz 

V 

/l  -  a2  -  (32  +a  +  p  Z_JlZ 


( 

du 


—  dxdy 

V  Ay 


=  J  U0(x,y, 0)  J  exp  [ikzF(a,fi)\  d 
x,y  a,P 


^ j  dxdy. 


(Bl) 


The  function  F( a,  (3)  is  expanded  in  a  Taylor  series  about  the  point  (a0,  (30)  at  which  its  first  derivatives  are 
zero.  Setting 


^Oc(aO’Po) 


(B2) 


and  similarly  for  (3,  leads  to 


and 


«0  = 


x'  —  x 


r 


Po  = 


/-y 


The  second  derivatives  of  F{ a,  (3)  at  (a0,  Po)  are 

Fact  («0  ’  Po )  =  — 


r[(x'-x)2  +z2] 


Vl-«o-Po  ^I-Oq-Po 


(B3) 


(B4) 


(B5) 


(B6) 
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^pp(aO’  Po)-- 


r\(y'  -  yf  +  z  ] 


(B7) 


and 


^ocp(aO’  Po)-  _ 


a0Po 


r(x'-x)(y'-y) 


1-ag  _Po 


(B8) 


The  following  quantities  will  be  needed  below: 


nao.  Po)  =  ^  +  =  r 

r  z  Z  Z 


(B9) 


and 


^aa (ao »  Po)%(«o.  Po)  -  ^p(a0’  Po )  — 


(BIO) 


The  Taylor  series  expansion  through  second  order  of  F(a,  (3)  about  the  point  (a0,  |30)  is 


F(a,(3)  =  F(a0,  (30)  +  —  Taa(a0,(30)  (a-a0)2  +  —  Tj3p(a0,p0)  ((3-(30)2 

+  ^ap (ao » Po )  (oc-a0)  (P-Po)- 


(B 11) 


The  a  and  (3  integrals  in  Eq.  (Bl)  can  now  be  evaluated.  Observe  that  Faa  =  - 1  Faa  I  ( Faa  is  always  nega¬ 
tive)  and,  for  ease  of  notation,  drop  the  (oiq,  |30)  argument  from  the  quantities  Faa,  Fpp,  and  Fa^.  The  first, 
second,  and  fourth  terms  in  Eq.  (Bll)  are  inserted  into  the  exponent  in  Eq.  (Bl),  and  the  integral  over  a  is 
evaluated  by  completing  the  square  in  the  exponent  (the  symbol  ±  in  Eq.  (B 12)  means  add  and  subtract,  not 
add  or  subtract): 


exp  [ikzF( a0,  (30 )]  J  exp  l  ikz 


=  exp  (ikr)  \  exp 


1  2 

~z  I  ^aa |  (®  —  ^ap (®  —  *-*0  )(P  —  Po ) 


ikz  ,  „ 
- Fr 


aa 


"a 

.2 


(a - a0 )2  - 2  r^&T(a-a0)(p-p0)±  — ^-(p-p0)2 


aa| 


[  aa 


=  exp  (ikr)  exp 


=  exp  (ikr)  exp 


^TT^(P-Po)2 


2  \F 

Z.  \1  f 


aa| 

2 


J  exp 

a 


ikz .  „ 
— —\Fr 


aal 


a-ao-nr^  (P-Po) 


-|2 


aa| 


\d 


V^T®-Po)2 

2  \Fn 


aal 


1  -i  1 


V2  ^ 


z\F, 


(B12) 
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The  integral  over  (3  in  Eq.  (B 1)  can  now  be  carried  out  by  adding  the  third  term  in  Eq.  (B 11)  to  the  exponent 
in  Eq.  (B 12),  and  again  using  Faa  =  - 1  Faa  I : 


.  exp  ( ikr )  J  exp 

vM^aal  P 


-t  ^(P-p0r+y  %(P-P0r 


1  -i  1 


^2 


( 

ikz 

Fap 

2  V 

E7 

aa| 

V 

J  exp 

/ 

.2  .  ikz 


aa 


1  -i  1 


^2 


exp  (ikr) 


oca 


p 

\-i 

VT 


f  r  rr  j-,2  \ 

^aa^pp  <;xp 


1  aa| 


(P-Po)2 


f 

1  aa 

^(^aa^pp 

-fap) 

IZ 

hi 


exp  (ikr). 


Eq.  (Bl)  can  now  be  written  in  the  desired  form: 

j  j {  /  /  \  *  r  TT  ,  exP  (ikr) 

U(x  ,y  ,z)  =  -  -  J  U0(x,y, 0) - - - 


X 


dxdy 


x,y 


=  -7  |  U0(x,y, 0) 

k  x,y 


r 

exp  (ikr) 


cos  %  dxdy . 


(P 

u 

FP 


X 


(B13) 


(B14) 


